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// for each pixel in the NxM image 

for(x=l;x<N;x+=l) 

{ 

for(y=l;y<M;y+=l) 
{ 

// Consider I(x,y), I(x~l,y), I(x,y-l) f ... neighborhood of (x,y) 

// If some ifl(x-l,y-j) fall out of image boundaries, replace them by 0 

// Do not consider reference pixels (I(x-i,y-j) where ij=0) 

// Apply (6) to compute e(x,y) 

sumJcI-O; 

for(i=0; i<a; 

f 

for(j=0;j<a;j+=l) 
f 

if(I(x-i,y-j)==out of bounds pixel) I(x-i f y-j)-0; 
if(i~~0 && j=-0) continue; 
sumjd =sum_kl +kj *I(x-i,y-j); 

} 

} 

sumjd -sumjd> > 1 0; 
e(x t y) =I(x f y)~sum_kI; 

} 
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int RTJPEG: :CompressRT(BYTE* input, BYTE* output, double ratio, double maxtime, 

int bpp=l, bool rgb=false) 

{ 

int sizejorig = Size(input); II original image size in bytes 
// Set max and min quality values 
int q0=5; II 5% quality 

int ql =95; II 95% quality 

// Find respective compression ratios. Higher ratio corresponds to lower quality 
double rO = size_orig/Compress(input,output,qO, int bpp=l, bool rgb=false); 
II highest ratio 

double rl = size_orig/Compress(input,output,ql, intbpp=l, bool rgb —false); 
II lowest ratio 

// Return if proposed ratio is smaller than the smallest possible 
if(ratio<rl) return JPEG: . Compress (input, output, ql); 
II Return if proposed ratio is larger than the largest possible 
if(ratio>rO) return Compress (input, output, qO); 
II Start timer 
clock J start, finish; 
start = clockQ; 

II Start iterative process to estimate the quality value 
double duration=0; double r = ratio; 

while (duration<maxtime) 

{ 

II Use linear quality estimate 

int q = qO+(r-rO)*(ql-qO)/(rl-rO); 

II Alternatively, a bi-section estimate can be used in the previous line: 
//intq = (q0+ql)/2; 

// Update estimated corresponding compression ratio value 

r = size_orig/Compress(input,output,q); 

II Update compression and quality boundaries 

if(ratio>r) 

{ 

rl=r; if(ql==q) break; II convergence 

} 

else 

{ 

rO=r; i f(qO==q) break; II convergence 

qO=q; 

} 

II Update timer 

finish = clockQ; duration = (double) (finish - start) / CLOCKSJ>ER_SEC; 

} II end of iterative process 
// Compress 

return Compress (input, output, q); 
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// It is assumed that ApplicationEntityList ael and DICOMObject d variables 
// have been already initialized based on the present network 
// configuration and the data to be sent. 

// Find current remote application entity, to which data must be sent 
ApplicationEntity ae = ael. GetCurrentEntityQ; 

II Create PDU connection 
PDUpdu; 

II Start network performance timer 
NetworkTimer nt; 

II Find data size to be sent 
nt GetDataSize(d); 
nt.StartTimerQ 

II Send data to current application entity 
pdu.Send(d,ae); 

II Stop the timer, evaluate observed network bandwidth, and store it into ae 

nl EstimateCurrentBandwidth(ae); 
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// It is assumed that ApplicationEntityList ael and DICOMObject d variables 
// have been already initialized based on the present network 
// configuration and the data to be sent. 

// Compress DICOM object to accommodate current network bandwidth 
int dsize — d. GetSizeQ; 

int ntime = ae.GetDownloadTimeQ; II maximum allowable download time 

int ctime = dsize / ae.GetCurrentBandwidthQ ; //time we'll spend with current bandwidth 

iffctime > ntime) II low bandwidth, we need compression 

{ 

iffae. CanAcceptLossyQ) II we can use lossy compression 

f 

// Compress with maximum ratio allowed 

double ratio = min(ctime/ntime, ae.GetMaximumCompressionRatio()); 
d. Compress (false, ratio); 

} 

else iffae. CanAcceptCompressedQ) II only lossless, try our best 

{ 

d.Compress(true, /*ignored*/L0); 

} 

else II no JPEG compression allowed, do nothing 

{ 

} 

} 

else II network is fast enough, do not have to compress 

o 

II Create PDU connection 
PDUpdu; 

II Send compressed data to current application entity 
pdu.Send(d,ae); 
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class AccurateTimer 

{ 

private : 

int Initialized; 

int64 Frequency; 

int64 Begin Time; 

public : AccurateTimerQ // constructor 

{ 

II get the frequency of the counter 

Initialized = QueryPerformanceFrequency( (LARGE _INTEGER 
*)&Frequency ); 

} 

void StartTimerQ II start timing 

{ 

iff ! Initialized ) return 0; II error - couldn't get frequency 
// get the starting counter value 

QueryPerformanceCounter( (LARGE JNTEGER *)&BeginTime ); 
return 

} 

double EndTimerQ II stop timing and get elapsed time in seconds 

{ 

if( ! Initialized ) return 0.0; II error - couldn't get frequency 
// get the ending counter value 
int64 endtime; 

QueryPerformanceCounter( (LARGE ^INTEGER *)&endtime ); 
II determine the elapsed counts 

int64 elapsed = endtime - BeginTime; 

II convert counts to time in seconds and return it 
return (double)elapsed / (double)Frequency; 

} 
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bool SoundRecorder: :Record(int sf, bool stereo, int max_time=300, int 

max_size= 1000000, int format=WAV) 

{ 

DWORD dwReturn; 

II Open a waveform-audio device with a new file for recording. 
if(!OpenMCI()) return; 

MCI RECORD FARMS mciRecordParms; 

II Set recording parameters here as fields in mciRecordParms 

mciRecordParms. dwFrom = 0; 

mciRecordParms. dwTo = maxjime; 

mciRecordParms.dwCallback = (DWORD) (this->GetSafeHwnd()); 
II Record 

if (dwReturn = mciSendCommand(m_Device.wDeviceID, MCI RECORD, 
MCI_FROM | MCIJVOTIFY, (DWORD) (LPVOID) &mciRecordParms)) 

{ 

AfxMessageBox(ErrorMCI(dwReturn, "Recording error: ")); 
mciSendCommand(m_Device.wDeviceID, MCI CLOSE, 0, NULL); 
return false; 

} 

return true; 

i 
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void SoundRecorder: :Insert(DICOMObject& dob) 

{ 

BYTE* sound; 

if(GetFormat() != MP3) II we have to encode the sound 
{ 

SoundEncoder se; 

if(se.Encode(GetSound(), GetFormatQ, sound) 

{ 

return false; II we failed to encode 
//Store sound bytes into odb object 



return true; 

l 
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Find Ri histogram: Remove outliers 
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